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ABSTRACT 

Analysis of five stable cases of the smoke pltmies that 
originated In eastern Cabo Frio (22^ 59*S; 42^ 02'U)» Brazil using 
LANDSAT Imagery Is presented for different months and years. From 
these Images the lateral standard deviation (o^) and the lateral eddy 
diffusion coefficient (K^) are obtained from the formula based on 
Taylor's tiMory of "diffusion by continuous moment". The rate of 
kinetic energy dissipation (c) is evaluated from the diffusion 
parameters and K^. Then, the vertical diffusion coefflclmt (K^) Is 
estimated using Uelnstock's formulation. These results agree well with 
the previous experimental values obtained over water surfaces by 
various workers. Values of c and stow the weaker mixing processes 
In the marine stable boundary layer. The data sample Is apparently too 
small to Include representative active turbulent regions because such 
regions are so Intermittent In time and In space. These results form a 
data base for use In the development and validation of mesoscale 
atmospheric diffusion models. 

1. Introduction 

There Is a strong tendency In the atmosphere and oceans 
for the development of layers In which the local gravitational 
stability of the air or water is either very large or very snail. With 
larger stable layers In the atmosphere or In the bod'ss of water there 
art often found sublayers, sometimes very ninnerous. In which neutrally 
stable layers are sandwiched between thin stable layers (Stommel and 
Fedorov. 1967). Some of this Is understood to some extent, but nany of 
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these layering phenonena still have no acceptable explanations. It Is 
evident that the dispersion of atmospheric pollutants In such stable 
layers Is limited by tte height of the Inversion lid. 

Considerable teterest has been evidenced In the 
properties of t1» Idealized planetary boundary layer contributing to 
the dispersion of trace naterlals Introduced into It. While the list 
of technologically significant atmospheric dispersion problems Is very 
long, those dealing with the fate of pollutants from Industrial 
processes have received special attention In the last decade as part 
of an effort to evaluate and predict air quality. Efforts to specify 
atmospl^rlc eddy diffusion coefficient for use In the models have 
centered f^lmarlly on an analysis of different studies (Tennekes, 
1973; Bauer. 1974; Mellor and Yamada. 1974; Yamada and Nellor. 1975; 
Zenan and Tennekes. 1977; Freeman. 1977; Bauer et a1.. 1978; Remsberg. 
1980). We mention some of the meteorological properties involved In 
this Important and difficult problem: wind speed, thermal 
stratification, wind shear and turning. These properties affect the 
dispersion of pollutants through the advective transport of the wind 
field and the diffusion which the pollutant experiences due to the 
turbulent motion of the atmosphere. 

In this paper we Investigate the vertical eddy diffusion 
coefficient In the atmospheric boundary layer from the LANOSAT inagery 
of the smoke plumes. It will be shown that che diffusion coefficient 
values derived under some approximations, with relatively simple 
expressions account qualitatively In a reasonable way for the above 
effects. 
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2. Theoretical considerations 


The relationship between the apparent lateral ed<1/ 
diffusion coefficient, (a* s*^), and the lateral standard 
Aviation, Oy (a), along the y- axis was derived by Batchelor (1949): 


•V 




(I) 


where t (s) Is the tiae. This expression, with Taylor's (1921) 
hypothesis becomes 




( 2 ) 


where u (a s*‘) Is the mean wind speed along the x>ax1s. The turbulent 
ener^ dissipation rate e (ro‘ s*’). Is estimated from the atmospheric 
turbulent diffusion study of Richardson (1926): 


Ky • 



(3) 


where here Is the scale of the diffusion process. The above 
treatment has already been presented In a detailed form elsewhere 
(Viswanadham and Torsani, 1982). 


Determination of the vertical standard deviation 02 ^*"^ 
and the corresponding vertical eddy diffusion coefficient, (m^ 
s“^), depend on measurements of concentration at different heights 
and hence are less reliable than those of o’ and Ky. Ozmidov (1965) 
and Lilly et a1. (1974) have proposed the following expressions for 
as 
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K, « 0.1 — ^ (Oznldov) (4) 

2 N* 

and 

K_ « — - — (Lilly et al.) (5) 

* 3N» 

in the ocean and in the stratosphere* respectively. Here* N (s*^) is 
the Brunt-Vaisala frequency. Eq. (5) for was based on the assumption 
that the flux Richardson's number (R^) is 1/4 and that the energy 
production is locally dissipated. The vertical turbulent diffusion 
coefficient in a stably stratified fluid is also derived analytically 
by Meinstock (1978)* i.e.* 

It » 0.81 . (6) 

'2 

Eqs. (4) and (5) are similar to the expression (6) derived by 
Weinstock (1978) except that Weinstock uses a coefficient of 0.81 
rather than 0.1 and 1/3. Since* theoretically most occurrences of 
turbulence should have small enough scales in the inertial subrange* 
it follows that Eq. (6) might be fairly generally applied to stable 
atmospheric conditions (Weinstock* 1978). Weinstock's derivation does 
not require that the flux Richardson number be known or specified. 
Recently* this expression was tested by experiments of Caldwell et al. 
(1980). Their measurements and analysis suggest that Eq. (6) is «1id 
for strong or intermediate stratification (actually* for Cox nunbers 
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less than 2S00) and Is not valid for Meak stratification (Cox nunbers 
exceeding 2500). fact that Eq. (6) cannot be valid for weak 
stratification Is evident from Its form, which shows a singularity at 
N - 0. 


Incidentally, Welnstock (1981) has also derived another 
expression for that applies to all stable stratifications (all 
positive N^). This derivation Is based on one Improvement of the 
derivation of Eq. (6). In the derivation of Eq. (6), It was assumed 
that the vertical velocity spectrum Is "cut off" at the buoyancy wave 
number regardless of how weak the stratification Is. However, a cutoff 
Is not valid when the stratification Is very weak (N 'x. o), as the 

stratification Is then simply too weak to Influence the spectrimu 
Hence the Improvenent Welnstock (1981) made Is to account for general 

stable stratification without using a cutoff at all. Eq. (6) gives a 

larger diffusion coefficient, by a factor of 2.4, than Eq. (5). 
Welnstock (1978) presented a detailed explanation to this fact. He 
suggested that the smaller ralue of obtained by Lilly et a1. might 
be attributed to the fact that turbulent fluxes were neglected (i.e, 
production was assumed to be local). This led to a smaller value of 
(■ 1/4) and a smaller \m1ue of K^. 

From the previous discussion we know that Eq. ((>) Is 

valid' for strong or Intermediate stratification. So, It Is better to 
compare the critical wave number k^ (m"^) of the energy-containing 
region of the spectrum with the microscale wave number km (m'^) of 
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Kolmogorov (1941). An expression for the critical wave number, k^, ms 
derived theoretically by Dougherty (i‘»61) as 

k . (7) 

” e dz 

where 6 (K) Is the potential temperature. Eq. (7) provides an 
estimate of the Mve nimber at which buoyancy cannot be neglected (see 

Welnstock, 1981; Fig. 1, p. 9927). The Kolmogorov's microscale wave 

nund>er, k|^. 1s 

V-/'. (8) 

Here, v (m^ s'^) Is the klnemtlc viscosity of air. Thus,the Isotropic 

Inertial range exists when kj^ > k^^ or Lj^ < L^^, where L (ro) Is the 
Mvelength (■ k'^). If, on the other hand, k^j < k^ or > Lq, there 

Is no Inertial range In the usual sense. In other words, we have weak 

stratification (I.e., most turbulence scales are not Influenced by 
buoyancy). In case, when k|^ < k^, we should not use Eq. (6) to 
calculate K^ In a stably stratified atmospheric flows. Assuming an 
Inertial range turbulence of the atmosphere, the values of K^ are 
obtained from Eq. (6) In the preser.^ study. 

3. Site and data sets 

Fabrics Naclonal de Alcalls (National Chemicals 

Industry) Is located on the east shore of Cabo Frio 120 km east of Rio 
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de Janeiro city In the state of Rio de Janeiro, Brazil. The Cabo Frio 
area (22° 59'S; 42° 02 *W) of eastern Brazil Is a flat coastal region 
bounded north'-west by land and the rest by Atlantic Ocean. The coast 
Is straight unconpllcated and typical of much of the northern and 
southern seaboard from Cabo Frio town (see Viswanadham and Torsani, 
1982). 


I 

\ 1 


:■) 
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The LANDSAT data were obtained during the period April | 

1975 to July 1978 at the above site. In all, 5 LANDSAT Images have 
been obtained during this period showing stable conditions. Figs. 1 
and 2 are typical examples of the LANDSAT Inages. The processing of 
the LANDSAT Images, and the determination of o and K from the 

Jr Jr 

LANDSAT Ima^s are presented elsewhere (Viswanadham and Torsani, 

1%2). The mesoscale meteorological conditions existing during the 
period (1975 - 1978} of interest could be approximated from available 
data of two agencies. Due to the lack of adequate meteorological 
measurements for Inland and oceanic areas of the Cabo Frio region at 
the time of LANDSAT imagery, the upper air observations data were 
taken frcxn the International airport of Galeao, Rio de Janeiro, Brazil 
(Fcrca /^rea Braslleira, 1974 - 1978). Surface meteorological 

parameters are available at a weather station very close to the 

Industry (Institute Naclonal de Meteorologla, 1968 - 1978). The 
radiosonde soundings for each LANDSAT Imagery were plotted for the 

three<klay period (one day before and one day after the passage of 

satellite) along with the surface data. Then, the Interpolated mean 
wind speed was obtained at 76-m stack height level. In Table 1, the 
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TABLE 1. Reference nurtber of LANDSAT Imge cases (n) and various Meteorological paraaeters. 
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base and top of the Inversion heights are presented along with other 
details. The LANOSAT Images have been characterized by a gradient 
Richardson number (RI) computed from radiosonde and surface 
observations. The RI values given In Table 1 are approximate due to 
the Interpolation Involved In the data sets. However, on the basis of 
RI values, we have five stable cases (Table 1). 

The lapse rates vfere estlmted over land <jhie to lack of 
marine meteorological observations. Iktfortunately, don't ha^ . 
evidence that the determined lapse rate over land Is a reasonable 
estimate for the lapse rate In the boundary layer over water. There 
are probably more Important differences In msoscale cllmte In the 
200 km east of the Cabo Frio region. ITie data are stmmarized In Table 
1 and do In fact suggest some meaningful Interrelationships (see 
Vlsvmnadiuim and Torsani, 1W2). The LANDSAT Image cases show that the 
smoke pitmies are very long and well (tefined except tte case n>8 (Table 
2). From these long and well defined visible pitmm characteristics of 
the LANDSAT Images, we nay have some evidence that the smoke plumes 
are associated with stable lapse rates the mrlne boundary layer. 

In the rwrmal lower atmosphere N* Is of the order of 10'** s'*. 
OeMrally, does not show abrupt variations in Its order of 
magnitude from near Inland through amrlne atmosphere. Therefore, In 
the present study the estimated lapse rate over land does not create 
abnormal differences In the calculations of vertical diffu, i';n 
coefficient In the marine atmosphere. 


Estimation of the rate of kinetic energy dissipation (e) and vertical eddy diffusion 
coefficient (K ) from smoke plumes of the LANOSAT Imagery. 
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4. Results 

As suggested In the Introduction, the phenomenon of plume 
spread can be quite complex unless some simplifications can be nade by 
careful scrutiny of the physics of the problem. Therefore, In this 
section Important deductions are made from the LANDSAT Imagery of smoke 
p1ume» and established physical laws of fluid n^chanlcs. On the basis 
of RI values, we classified the LANDSAT Inages Into 5 stable cases 
(I.e, n s 5, 7, 8, 10 and 11, In Table 1). Then the measurements of 
plume geometry provided the estimate of lateral standard deviation 
and the lateral eddy diffusion coefficient for 5 cases (Viswanadham 
and Torsani, 1982). 

Monin (1959) suggested that the diffusion along the 
vertical In the lower layer of the atmosphere Is mostly due to 
small-scale turbulence. And, the Richardson law K (o^) a 
explained by Obukhov (1941) as a consequence of the similarity 
hypotheses of Kolmogorov (1941) for turbulence with a large Reynolds 
number. According to Kolmogorov's hypothesis, the rate of kinetic 
energy dissipation e Is the only parameter determining the turbulent 
regime within the Inertial range. If the diffusion occurs as a result 
of turbulent motions of scales within the Inertial range, then Eq. (3) 
Is obtained. In other words we obtain Richardson's law. Therefore, the 
parameter e by viscous process can be regarded as a good Indicator of 
the strength of turbulence. With the aid of theoretical considerations 
In the section 2, the LANDSAT data sets were used to make estlnates of 




- 13 - 


2*25!S ^ It 

^ on*- -vy 


c and the vertical eddy diffusion coefficient In 5 stable cases. 
The critical vAvelength and the Kolmogorov's microscale imvelength 
were also evaluated from Eqs. (7) and (8)» respectively. In Eq. 
(8), the value of 0.15 x 10"" m* s”‘ Is used for the k^ematlc 
viscosity of air v (Sutton, 1953, p. 38). The results are presented 
In Table 2 together with other parameters. In Table 2, mlues of e 
and show the weaker mixing processes In the mrlne boundary layer. 
The mean c values show how sensitive the thermal Inertial boundary 
layer can be to changes In wind. Insolation, etc. They must be 
regarded as tentative, due to the subjective nature of the 

Identification of turbulence categories, and the diverse oata sources 
used In ccwnputations. In addition. It Is of great interest to 

consider the nwnerical values of and Ln In Table 2. In all 5 

cases, Ljj Is much greater than L^j. Thus, the Isotropic Inertial 

subrange spectrum exists In all cases. This argument shows that the 
utllltatlon of Eq. (6) to calculate In these 5 cases Is well 
Justified. 


Although one might expect some \«lues for the parameter 
Kj, a check with actual measurements is more important. In Table 3, 
we have sixmerized results of from four available publications. 
Note that the reported MLlues of K^ In Table 3 are directly 
comparable with our results In Table 2. Lilly et al. (1974) have 
showed that e is subjected to large variations with underlying 
surface characteristics, height, latitude and season. Also, the 
measurements over water are difficult; because the spectra extend to 
very small scales. 
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TABLE 3 OUAimr 

Values of Kj from tlw llttraLure 


Source 

Location 

«• «■* 

Boll (196S) 

Over water 

0.05 - 3.50 

Bowden et al. (1974) 

Sea areas 



Irish sea 

10*» - i0*» 

Lyons (1975) 

Over water 

1.50 • 3.00 

Bond and Bryan (1976) 

Ocean area 

0.01 

Bresent results 

Ocean area 

0.01 - 3.42 


There are several possible sources of error in the 
above computations, ^ceptable errors vrauld have been encountered in 
all Instances In any case because of satellite |K>sition in relation 
to the study area* and the size of an individual grid cell compared 
to the magnitude of possible error (Ernest and Lyons. *974; Ernest, 
197S). Some error may arise from assuming an inertial range spectnmi 
to compute dissipation rate. Also, the vertical, lateral and 
longitudinal turbulence differed substantially fr«n the assumption of 
isotropic turbulence. There arc no direct atmospheric data against 
which to test these preditions. 


Over water, while the surface roughness tends to 
be- more nearly uniform than one could find at land sites, sharp 
variations in water surface temperatures, especially in near shore 
region, can mke the Interpretation of data exceedingly difficult. 



IS 


Furthtnwrt, during ptrlods of offthort f1^» IntMiso air mss 
modlflcition occurs. Thus, ts sir flows fron Und to cold ustor, not 
only <^s tht boundary ta^raturt change, but tha value of the eddy 
conductivity begins to decrease to a lower value. There Is a 
(testructlon of the MChanlsm Involved In creating high values over 
land, since much of the vertical eddy transport probably results from 
penetrative convection acting over and above the effect of mechanical 
turbulence. It Is this convective component that disappears very 
rapidly with fetch. The drop In from around 10 to 10"^ m* $*' (or 
less) probably occurs within the first few kilometers. The air's 
vertical diffusion rates should Indeed be low In such an environment. 
There Is little or no quantitative Information on fetch, wind 
duration, sea state, etc., with which we can categorize Oy In the 
present analysis. Their effects (If any) remain obscure In the data 
scatter of Individual data sets (Tables 1, 2). In summary, the 
present data sets give evidence that the average vertical diffusion 
In the lower troposphere Is somewhat lower In the stable marine 
boundary layer. On the other hand, the vertical diffusion 
coefficients of smaller magnitude compared to ours are found 
appropriate when the quasi-horizontal processes are explicitly 
represented In the global atmosphere. 

6 Concluding reimirks 

NHSurements In the atmosphere even under the most 
favourable conditions could Nrdly be Interpreted directly In any 



particular case, even If e and the distribution of concentration could 
be measured accurately; because the turbulence Is certainly not 
sufficiently uniform over Its depth to permit uncritical application 
of a theory which depends on uniformity as well as local Isotropy. The 
reason why Kolmogorov's analysis seems to apply so well to 
Richardson's ^nthesls of all observa.t1on of diffusion has been 
touched on by Batchelor and others. In spite of the Importance of 
those diffusion coefficients for the quantitative treatment of 
atmospheric exchange processes. It must be stated that no theoretical 
approach attempting to Interpret these highly variable factors by 
characteristic quantities of the turbulent and convective motions have 
become known up to now. Hence, we have to be satisfied with empirical 
data derived from these suitable measurements. 

Finally. It Is obvious that this study Is a first 
attempt at a statistical description of the turbulent microstructure 
In the marine boundary layer rather than a definite study. Until the 
nature of the microstructure is proven conclusively our understanding 
of the mechanisms producing the mixing events will be limited. 
Therefore, much more needs to be known about the physical processes of 
stratified turbulence before its signature can be properly identified 
In the marine boundary layer. In the imrine atmosphere, however. 
Ingeneous observational techniques had to be applied before any 
rele\^nt results could be obtained. 
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Fig. 1 - The June 23, 1976 LANOSAT Imge at 1130 GMT of smoke plume of 

the Cabo Frio region drifting NE iMrd across the Atlantic 

Ocean. The air was so stable that tte plume could be seen 

extending for 130 km with u ■ 4.1 m s~^ and ~ > 5.6 K km*^. 

3z 

But, over water fetch the air temperature mss slightly colder 
than the ocean surface tei^rature. The fomatlon of cumulus 
cloud streets are visible over the ocean In right side of 
photo. Image photographically enhanced from original 35 m m 
slide. 

Fig. 2 > As In Fig. 1 except for 1205 GMT July 1, 1978. The air was so 
stable that the plume could be seen extending for 130 km with 
u ■ 4.6 m s’^ and ~ « 2.7 K kra"^. A strong clockwise surface 
water current is visible on right side of the Island. 


